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We study Bianchi I type brane cosmologies with scalar matter self-interacting through combina- 
tions of exponential potentials. Such models correspond in some cases to inflationary universes. We 
discuss in detail the conditions for accelerated expansion to occur: in particular, we show that the 
condition which is necessary and sufficient for inflation in the relativistic version of the models is not 
sufficient in the brane case. Another peculiar feature of the models is that the relationship between 
the value of the scale factor at the beginning of inflation and the equation of state is very different 
from what one finds in the relativistic framework. We also analyze the influence of the value of the 
anisotropy and the brane tension, and show that the associated effects become negligible in the late 
time limit, those related to the anisotropy disappearing earlier. This study focuses mainly on single 
field models, but we also consider a generalization yielding models with multiple non-interacting 
fields and examine their features briefly. We conclude that, in the brane scenario, an increase in the 
number of fields assists inflation, as happens in general relativity. 
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I. INTRODUCTION 

According to the most promising candidates for a 
quantum gravity theory, we are hving in more than four 
dimensions , and gravity is effectively four dimensional 
only at low enough energies. This has inspired scenar- 
ios in cosmology different from the standard one. The 
"brane world" picture is the most recent of them all 0, 
and assumes that ordinary matter is confined to a hy- 
persurface with three spatial dimensions embedded in a 
multidimensional space-time, allowing gravity act in the 
fifth dimension. 

The "brane-world" proposal has motivated a revision 
of the predictions to be tested by future observations. 
Even if in a first approach isotropy is assumed, one will 
immediately find Q that the evolution of the Universe 
in this new scenario is rather peculiar, at least at early 
times, because of the modifications arising in the Hubble 
equation. Given this, one may wonder for instance what 
aspects of cosmological inflation should be changed to fit 
in this alternative description. Clearly, one of the most 
important problems to address within this framework is 
that of the initial conditions. Recently, it has been sug- 
gested that anisotropic spacetimes typically expand 
faster than isotropic ones. This study was grounded 
on a qualitative analysis of the evolution equations for 
the scalar field (Klein-Gordon equation) and the average 
expansion rate (generalized Raychaudhuri equation) for 
Bianchi I type cosmologies. 

In this spirit, we have investigated the same equation 
set, together with its solutions, just assuming at the be- 
ginning that the potential depends on the scale factor 
in such a way that it is compatible with inflation. In 
particular, we have found that, under some assumptions, 
the matter source can be cast in the form of a scalar 
field self-interacting through a combination of exponen- 
tial potentials. Such models were first considered in the 



context of inflation nearly two decades ago p| , and have 
been paid much attention along the years because of sev- 
eral reasons. First of all, there is solid motivation that 
makes them physically appealing, for they appear in four- 
dimensional effective Kaluza-Klein theories arising from 
compactification of higher dimensional supergravity or 
superstring theories (see 6] for a review). In addition, it 
was shown long ago that FRW power law models, which 
are driven by exponential potentials, behave as attractors 
in the sense that models with broad disparity in their ini- 
tial conditions become power law solutions at late times 
0, Last but not least, exponential potentials have 
proved very useful in general relativity for providing new 
exact solutions because they add a rather small degree 
of non-linearity to the field equations 0] . Unfortunately, 
the equations one has to solve in brane cosmology are 
intrinsically so non-linear that, even if one considers just 
exponential potentials, the task of finding exact solutions 
is a hard one, except in very specific cases. 

With the aim of getting round this difficulty, and as 
mentioned before, we will not assume from the start that 
we have an exponential potential or a combination of 
such terms, but rather end up having it after making a 
certain ansatz regarding time evolution. 

As expected, if the expansion is fast enough then 
anisotropy exerts less influence on the expansion than 
the quadratic energy-momentum corrections due to the 
extrinsic curvature. It follows then that anisotropy does 
not prevent inflation, as one also gets in the framework 
of general relativity [wL Hit . In the late time regime, 
the effect of anisotropy and extra dimensions are negli- 
gible, and the solutions behave as if they were isotropic 
and purely relativistic. 

The qualitative behavior of Bianchi I type brane cos- 
mologies has been studied using dynamical systems tech- 
niques by several authors. Barotropic perfect fluids were 
considered by Campos and Sopuerta in the absence [l^ 



2 



and presence [Tj| of bulk effects; whereas Goheer and 
Dunsby [l3 considered models with a single scalar field 
(see also for an earlier study along the same lines 
focusing on isotropic cosmologies). According to those 
studies, late time isotropization is an inherent feature to 
the models. This idea of inflation acting as an isotropiz- 
ing agent, that was put forward years ago by Wald 1 171 . 
has been recently revisited in the brane scenario [ig 
(see also 0, |23| for a discussion on the isotropization 
of anisotropic relativistic cosmologies with an exponen- 
tial potential). Nevertheless, we have recently shown 
that the brane might interact with the bulk in a way not 
compatible with late time isotropization. Here, however, 
we will discard the bulk effects and Wald's argument will 
still apply. 

In another recent study Bianchi I type exact so- 
lutions in parametric form were found for barotropic 
perfect fluids and several scalar matter models. It was 
shown that for exponential potentials the simple relation 
(j) oc logt will not hold at all times. This does not ex- 
clude, though, the possibility of having this behavior at 
some epoch, and that is indeed the case for the late time 
regime of our solutions. We also show here that the ef- 
fects of the extra dimensions make the models isotropize 
faster, this being linked to an increase in the Hubble rate. 

Another interesting aspect of the problem is that, typi- 
cally, exponential potentials only render inflation if their 
slope is small enough. Since our solutions mimic their 
relativistic counterparts late in their history, it is not 
surprising that the same restrictions on the values of the 
slope will hold for our models too. This is discussed 
in Section 3. However, one could still have inflation 
with sloppy potentials if there are enough non-interacting 
fields. This is the so called assisted inflation proposal. 
We also analyze this possibility in the last section, and 
find that inflation is likelier in such multifield models be- 
cause the Hubble factor grows as the number of fields 
increases. 

The plan of the paper is as follows: in Sec. 2 we study 
the equations of motion for the scalar field and the gravi- 
tational field equations discarding bulk effects and under 
the assumption of a Bianchi type I geometry. In Sec. 

3 we consider solutions such that their potential energy 
decreases as the models expand, and find approximate so- 
lutions with power law inflation at very late times. Sec. 

4 is devoted to discussing in detail the inflationary be- 
havior, paying special attention to the role of anisotropy 
and the corrections from extra dimensions. In this sec- 
tion, we show that the value of the expansion factor at the 
beginning of inflation depends on the equation of state 
in an peculiar way, if we compare the results with what 
one gets in the relativistic. This suggests that the evo- 
lution of brane cosmologies is in some aspects very dif- 
ferent from what is regarded as standard, and that the 
subject deserves further investigation. In Sec. 5 we turn 
to consider multi-field generalizations and examine their 
features. We show that it is possible to model assisted 
inflation in the brane setup too. Finally, in Sec. 6 we 



draw our main conclusions, and review future prospects. 



II. GEOMETRY ON THE BRANE AND FIELD 
EQUATIONS 

We assume that the 5-dimensional line element giving 
the geometry of the bulk is of the form 



ds^ — gAsdx dx 



(1) 



where A,B = 0,1,2,3,4. Tildes wiU be used to de- 
note the bulk counterparts of standard general relativis- 
tic quantities. If we now denote by the unit vector 
normal to the brane, then the metric induced on that 
hypersurface is qab = Qab — riAns- It is convenient to 
define a new coordinate % such that the brane is located 
at X = and UAdx^ = dx- This renders the line element 
in the alternative form 



ds^ — dx^ + gfj,i,dx^dx'^ 



(2) 



with fi,!/ — 0,1,2,3. The five-dimensional bulk field 
equations are 

Gab = -~K^[^~9AB + 6{x)iX9AB ~ Tab)\- (3) 

Effectively, the Einstein equations on the brane take 
the form 



where 



TT 

12 



E„ 



(4) 



rji Ctrp 



24 



(5) 

A = k2(A + k2AV6)/2, k2 = K^A/e, T = T^, and 
Eab = CAiBjn^n'^ ■ Here Caibj is the five dimensional 
Weyl tensor in the bulk, whereas T^^ and A stand for 
the values on the brane of the energy-momentum tensor 
and the vacuum energy respectively. The equations were 
obtained under the assumption of Z2 symmetry, and con- 
finement of the matter fields to the brane, in agreement 
with the brane scenario devised by Horava and Witten 

It has been shown that i?^^ can be decomposed like 



(6) 

where is the four velocity of an observer on the brane, 
U is the effective non local energy density on the brane, 
Vfiu is the effective nonlocal anisotropic stress, and 
is an effective non local energy flux on the brane, which 
vanishes identically for a Bianchi type I brane. In addi- 
tion V^u and are subject to the constraints 



(7) 
(8) 
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Let us consider now four-dimensional Bianchi I type 
spacetimes. The line element is given by 



a\[t)dx^ + al(t)dy^ + al{t)dz^ 



(9) 



As usual, we define the expansion rates along the three 
spatial directions as Hi = hi/ ai for i = 1, 2, 3. The shear 
tensor will be denoted as and 



a^, = V(i?. - H f, 



(10) 



where 3H = Hi + H2 + H3. For this metric, Eq. Q is 
identically satisfied. Moreover, since there is no evolution 
equation for Vfj,i, on the brane, one usually assumes that 
it is either null or satisfies the less restrictive condition 
that (J>""P^,y = (It has been proved that the inte- 
grability conditions for = 0, V^v = imply spatial 
homogeneity, so this is consistent with having a Bianchi 
I metric on the brane.). Both assumptions on Vfj^i, trans- 
form Eq. ((SJ into a very simple evolution equation for U, 
namely 



U + AHU = 0. 



(11) 



The solution of the latter corresponds to non local en- 
ergy evolving like radiation, i.e. lA = Uo/a^, with Uq a 
constant that can be null. As we see, the joint choice of 
V^ij — and U — is consistent on the brane. Whether 
this is also fully consistent for the bulk remains an open 
question, but we follow other authors in making these 
arguably strong assumptions 0, . This allows gaining 
considerable insight in the changes on the evolution of 
cosmological models associated with the modification of 
the Friedmann equation, which might be otherwise con- 
cealed by bulk effects. 

Under the assumptions made, the gravitational field 
equations and the Klein-Gordon equation read 



HiH2+HiH3 + H2H3 = A + >i^p+—p^ (12) 



Hl + H2H^ + Hl + H2 + H^ = k-n'p-^{2p + p) k\ 
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(13) 



Hl+HiHi + Hl + Hi+Hi = K~t^p-^{2p + p) k\ 

(14) 
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Hf + HiH2 + Hi+Hi+H2 = A-K^p--^{2p + p) k\ 
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• dV 

3H(j>+— = 0, 



(15) 
(16) 



according to the customary interpretation of irrotational 
scalar fields in terms of a perfect fluid. 

Due to the isotropy of the matter source, one can ob- 
tain constraint equations which provide a relationship ex- 
pressing the expansion rate along one direction in terms 
of the other two and their derivatives. 

Even though the solutions to the system H12(l - (|15|l in- 
clude locally rotationally symmetric models (LRS), i.e., 
models with the same expansion rate along two direc- 
tions, we will not consider them any further. Thus, 
in what follows we will assume Hi ^ H2, H2 ^ H^, 
i?3 ^ Hi, so that the models admit only three isome- 
trics. The constraints arising from having identical pres- 
sure along the three spatial directions are 



Hi + H2 + -ffs 



H2 — Hi 
H2 — Hi 



H^ — H2 
H3 — H2 



Hi — H3 



Hi — Hz 
(19) 

This allows to distinguish two different cases. The first 
one corresponds to H2 ^ i?3. Hence Hi ^ H2 and 
H^ 7^ Hi. Integration of Eq. H19|) yields the following 
relationships: 

Hi - H, = ^{Hi - H2) = ^^{H2 - H3), (20) 

where ^ 1 is an arbitrary integration constant. The 
second case arises when the Hubble factors along any 
two directions differ only by a constant, then necessarily 
Hi + H2 + H3 — and we get a model which does not 
evolve on average. For simplicity we will discard these 
cases too. 

It is convenient to write the field equations using the 
following definitions: 

P = Hi-H2, (21) 
Q = Hi-H:i. (22) 

We will also use the following constant: 

^ >4- (23) 

Note that, since we are not considering LRS models, we 
will have w ^ 1/3. Equations ((T^ - l(TC|) become 



3P 

A + n'^p+—p 



■ dV 
3H(f) + — = 0. 



(24) 
(25) 

2 - 3H^ - mP^, (26) 
(27) 



III. 



where we have set 



COSMOLOGIES WITH EXPONENTIAL 
POTENTIALS 



P=-^^ + V{cj,), 



(17) 
(18) 



In this section we look for solutions by means of an 
unconventional technique that consists in taking the av- 
erage scale factor a = (010203)^^^ as the independent 



4 



variable, instead of the time variable t. This procedure 
has been successfully used in related problems. 

Taking advantage of the insight gained in Ref. jl^l , we 
write 



(28) 



and make the change of variables dt — a^drj in Eq. (|27|l . 
so that 



:dV 



0, 



which in turn gives us the first integral 



1 



</)2 + y(0) 



6 



(29) 



(30) 



where c is an arbitrary integration constant. This is 
equivalent to saying that 



. ^ 6 /• , c 
PW = -fi / -aa+ -fi- 



(31) 



Combining H = d/a with Eq. (|25|l we get P — c'a^^, 
where c' is another arbitrary integration constant that we 
will set to one for simplicity. Now, using Eqs. H25|l and 
()26|l we can formally reduce the problem to quadratures, 
namely 



^/Sda 



(32) 



A + K^p{a) + —p{af + 



w 
Tie 



and 



V6 



p{a) - ^da 



A + K^p{a) + —pia)^ + ^ 



(33) 



12' 



where 0o are integration constants. We will choose 
the plus sign in Eq. (|32|l so that we have expansion as 
time grows, and will also make the same sign choice in 
Eq. (|33|l for reasons that will become clear in the lines 
below. 

Our main objective at this stage is to find expressions 
for F{a) so that they are suitable for depicting a given 
kind of inflationary behavior, and in particular that as- 
sociated with exponential potentials. We assume then 
that the expansion is driven by a scalar field striving to 
make its potential energy as small as possible. If the po- 
tentials have positive slopes, the only possibility is that 
the scalar field grows with time, and in particular with 
the scale factor. Therefore, we conclude that given our 
purposes we have to choose the plus sign in Eq. (|33|) . 

One simple way of satisfying the requirements we just 
made on the potential is to set 



F{a) = fa- 



(34) 



with n < 6 and / > 0. This gives 
V{a) = /a"-6. 



(35) 



Note that, as pointed out in Ref. [23, the consequence 
of the usual assumption A < allows for an arbitrary 
value of A, and in particular it can be chosen so that it 
vanishes or takes some other wanted value . Let us set 
A = for the time being, together with c = 0. We then 
get 



d<t) 
da 



= ±- 



V6 



(36) 



6/a" 



2na4- 



The relation V = Voe^^"^ will only hold if </) cx log a, so 
the question arises of whether that condition is compat- 
ible with Eq. (|33() . In order to find the answer we power 
expand d4>{a)/da. around a — 00 and keep only the first 
term in the series. It can be noticed immediately that 
the required condition <j) cx logo will only be satisfied 
asymptotically for large a if < n < 6. Therefore, the 
study will be restricted to just those cases. This would 
then give 



V^fe 



-n K(f) 



with 



■ logo. 



(37) 



(38) 



There are three exactly integrable cases, namely n = 
2,3,4, and the corresponding potentials turn out to be 
functions of e^'^^^^""^. For other values of n one can 
only give approximated expressions for late time lim- 
its. The leading terms in the solution 0(a) will be 
the -nj-dependent or the K-dependent ones depending on 
whether 0<n<3or3<n<6. 

With the aim of giving the expressions in a form as 
simple as possible we define the coefficients 



(39) 



24/k2 



and 



8nK2 



4riA' 



(40) 



where the letters s and h indicate that they are related 
to the shear and the brane tension respectively. We out- 
line the expressions for the scalar field and potential in 
Table ini 



n = 2 



n — 3 



n = A 



<n< 3 



3 <n< 6 



<^(a) 



1 + 2s2 + ^/a2 (a2 + 4s2) + 462 
-log ^ 



log ^ 



\/3 K 



1 + 264 + Ja^ (a2 + 464) +4s4 
— log 



S4 



2 ^6~^ , a"/^ + Va" + 4 s„ 
log 



fe 



— n /s:(}!) 



0/ \/6 — n 



2{6 — n)/n 



2 , a^-"^^ + Va«-" + 46„ 
log 



-\/6 — n 



/ 



TABLE I: Scalar fields and potentials for the integrable cases. The expressions in the last two rows are only valid in the large 
a approximation. 



We have chosen the integration constant (/iq in such a 
way that in the large a and large 4> limits we consistently 
recover the same expressions as in Eqs. H37|l and (I38II . 
Solutions similar to these were obtained by Barrow and 
Saich |2§| in a study of Friedmann universes containing a 
perfect fluid and a massive scalar field under the require- 
ment that the kinetic and potential energies of the scalar 
field be proportional. 

Before we go further, a remark regarding the choice of 
the value of A is in order. We have found three exact 
solutions for the A = case (see Table 1), but we can, 
by the same token, find another three for any non null 
value of A. This requires a little generalization that lies 
in including in the potential a term acting like a cosmo- 
logical constant. To this end we set F{a) — fa^ + ha^ 
which will obviously yield V{a) = fa"^^^ + h. Now, for 
the choice 



/i = -2 



3k2 ± V9k4 „ 3A;^4 



IV. KINEMATICS AND INFLATIONARY 
BEHAVIOR 

The subject in hands now is the kinematics of our cos- 
mological models. Let us first see if the energetic features 
of our cosmologies are suitable for inflation. A necessary 
condition for the occurrence of accelerated expansion is 
that the potential part of the energy density dominates, 
and for this to be possible the potential will have to be 
flat enough. Let us denote 



and 



Pkv 



Ppot 



1 



V, 



so that p = pkin + Ppot ■ Speciflcally, we have 



(41) 



Pki'. 
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■ Ppot ■ 



(43) 



(44) 



(45) 



one gets 




"-2 2na4-'^ 



(42) 



This can be obtained from Eq. H33I) upon the replacement 
— !■ + hn^ and the same applies for the expression 
giving dt/da. Therefore, the results for the previous set 
of solutions can be easily extended to any A 7^ value. 



This means that the kinetic and potential energy decrease 
at the same pace at late times (large a) if c 7^ 0, and at 
any time if c = 0. Moreover, if the latter holds we will 
have a barotropic equation of state p = {\ — n/3)p. This 
means we have found the exact potentials for radiation 
{n — 2 , p = p/3), dust {n — 3, p = 0), and cosmic strings 
{n — 4, p — — p/3) models. Note that, in the c ^ case 
the energy density could be dominated by its kinetic part 
at early times. 

It has to be pointed out that there is one more inte- 
grable case of Eq. H33|l . the stiff fluid model (/ = 0), 
which was found in Ref. 1231. 
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Typically, the inflationary behavior depends on the 
slope of the potential, but comparisons between po- 
tentials with different slopes must be done for identi- 
cal values of the scalar field. In our models we have 
d^V/{dwd(t)) > and d^V/idkd<j)) > 0. This can 
be checked using dV/d<p — [dV/da){da/d(l)) and Eqs. 
(|35|l and This would seem to indicate that the 

anisotropy and the quadratic corrections would make in- 
flation less likely. 

Nevertheless, the definitive answer is provided by the 
deceleration factor q, which is defined through 



and can be evaluated on calculating t — t[a) from Eq. 
(|32|l . In Table 2 we summarize the behavior of q at late 
and early times. 





t-to 


Q 


a ^ 


< n < 3 


/(6 — nja"^ 


5 — n 


n = 3 




2 


3 < n < 6 




2 


a ^ 00 


< n < 6 


(6 — n) \fj K 





TABLE II: Asymptotic expressions for t as a function of a, 
together with the corresponding deceleration factors. 

For large enough a we get <? = (4 — n)/2, so inflation 
will only occur at late times if and only if n > 4. This 
result is completely equivalent to that obtained in in 
the framework of general relativity. There is a persistent 
non inflationary behavior in the asymptotic regime of the 
models if the potential is too steep. For a ^ 0, however, 
inflation will by no means occur. 

Let us push a bit further the fact that the sign changes 
in q mark the transition between an inflationary and a 
non inflationary epoch. As follows from Eq. H46(l accel- 
erated expansion will occur if a > 0, which in our models 
holds if 

3 (n - 4) na^+'^fK^ + 3 (ti - 5) a^'^fk'^ - 2n^wa^ > 0. 

(47) 

We see that n > 4 is a necessary, but not sufficient, condi- 
tion for accelerated expansion, and also a sufficient con- 
dition for a ^ cxD, as proved above. Nevertheless, Eq. 



(|47|l has in store some interesting information that will 
highlight how peculiar inflation is in brane cosmology. 

Let us make stand for the value of the expansion 
rate above which the model expands. We will now study, 
by numerical techniques, how a-mi depends on the param- 
eters of the model. First, we will focus on the dependence 
on the equation of state. After that, we will turn to con- 
sider the effect of the magnitude of the anisotropy and 
the brane tension. 

Numerical calculations indicate that, in the relativis- 
tic limit, as is well known, the larger n, the smaller Oinf 
(i.e. the earlier inflation begins), whatever the value of 
Tu. Inspection suggests, however, that this monotonic 
behaviour may not hold for the brane case, because the 
term in Eq. H47|l associated with the extra dimensional 
corrections changes sign at rt = 5. Moreover, the values 
of Oinf in the vicinity of rt = 5 must be very much alike 
the values corresponding to the relativistic limit (because 
the extra dimensional corrections will be negligible) . Nu- 
meric evaluations show that for n close to 4 the behaviour 
of Oinf is almost identical in the brane and the relativis- 
tic case. However, as n = 5 is approached, the mono- 
tonic decrease ceases and a-mi begins to grow with n, till 
it reaches a maximum value before n = 5. Note that 
the value of the maximum of ainf does not depend on 
Tu. From that point on, there is again monotonic de- 
crease, but is more pronounced that in the relativistic 
case. Summarizing, the conditions for inflation are rather 
different in the brane and relativistic setups. Neverthe- 
less, several restrictions on the parameters of the mod- 
els have to be considered. On the one hand the poten- 
tial must remain much smaller than the fourth power of 
the four-dimensional Planck's mass Mp so that classical 
physics is valid. On the other hand, the five-dimensional 
Planck's mass Mp is typically much smaller than its four- 
dimensional counterpart. In Fig. 1 we have represented 
logOinf as a function of n, for several values of tu, within 
the ranges 4 < n < 5 and 1/2 < < 8 respectively, 
choosing Mp considerably smaller than Mp. 

We have seen, also by numeric means, that Oinf in- 
creases with the shear and decreases with the brane ten- 
sion^. However, those very changes in the parameters 
which act to delay inflation, make the Hubble factor grow 
in turn , and inflation can be sustained by potentials 
that would not be able to do so in a relativistic setup 
|30l |. Thus, both the anisotropics and the extra dimen- 
sional corrections seem to act locally in a way conducive 
to inflation. In order to determine whether this is also 
the case in a global way, it would be necessary to evaluate 
the number of e-foldings between the beginning and the 
end of inflation. So, we would have to flnd physical ar- 
guments allowing to know when inflation ends, and then 
calculate the amount of inflation achieved while acceler- 



The shear increases with zu, but the brane tension decreases as 
k increases. 
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FIG. 1: Plot of logai„f as a function of n for m = 0.25,0.5, 1,2,4,8, a(0) = l, / = 10"^°, n = Svr, and k = 8n x 10^. for the 
figure on the left and k = for the figure on the right. Here the Planck mass on the brane Mp has been set to one. Lower 
curves in each plot correspond to lower values of zu. 



ated expansion persists. 

Unlike what happens with the time at which inflation 
ends, we can determine for our models the time at which 
inflation begins. We have seen that the larger w or ii, the 
larger Cinf • But because of the way the Hubble factor 
behaves, the value of At associated with a given value 
of a grows too. Due to the lack of exact expressions for 
the t{a) curves for the inflationary cases (n > 4) one 
can only find the corresponding value of time by using 
approximate solutions or numerical integration; we will 
choose the second method. In Fig. 2 we have plotted the 
t{a) curves for of n = 4.5 and for the same values of Mp 
and zu as in Fig. 1; we have also indicated in the picture 
the value of Oint as a function of time. 



4-10 
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Plot of the t{a) curves for zu = 0.25,0.5,1,2,4,5 

-30 ^ r r _ in-30 



FIG. 2: 

a(0) = l, / = 10"^", n = 4.5, / = 10"^", k = Svr and k 
8tt X 10^. Here the Planck mass on the brane Mp has been 
set to one. Lower curves correspond to higher values of -oo. 
Each dot marks the beginning of infiation 

For the sake of completeness and comparison, let us 
investigate whether our models isotropize or not. We 
have seen above that the infiuence of the anisotropics in 
the inflationary behavior becomes less and less important 



as time goes by. The anisotropic shear parameter, which 
is defined as 



hear 



6H^ 



(48) 



refiects that too. Using Eqs. H21|l - H23|) . and after a little 
algebra, we get 



^shcar 



and according to Eq. H33() one finally has 



n. 



shear 



(49) 



(50) 



In the inflationary cases, since n > 4, we see that VLshear 
decreases monotonically with time, tending to a null 
value. Moreover, the larger n and k the more rapidly it 
decreases. Thus, although we have genuinely anisotropic 
models at early times, anisotropy gets completely dissi- 
pated in the course of the evolution. So, our models do 
isotropize at late times, and both faster inflation and ex- 
tra dimensional effects make the isotropization process 
more efficient. 



V. GENERALIZATION TO MULTI-FIELD 
MODELS 

Let us now take advantage of the results in the previ- 
ous sections to discuss how the situation would change 
if we had multiple fields. Assisted inflation jsj, in its 
simplest form, is realized with configurations of identical 
multiple fields with no interaction among themselves, but 
which self-interact through the same potential (see also 
33, 34] ) . Aspects of this proposal in connection with the 
brane scenario have been discussed in 133, IsM 123 . 
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If wc now make (j) denote each scalar field and Vi^cj)) 
each individual potential, we can obtain the multifield 
version of our models just by making the replacement 

K ^/Nk, (51) 
K ^ VNk, (52) 

where N denotes the number of fields. It can be imme- 
diately seen than the larger N, the larger the Hubble 
factor, so that the assisted inflation proposal is valid in 
the brane scenario too. 

The asymptotic individual potentials will become 

V ~ e-''V^(^-")'^. (53) 

Thus, an increase in the number of fields makes the po- 
tentials turn flatter, which means that late time inflation 
becomes likelier. 

The condition for having q < reads now 

3 (n - 4) n a''+"/ (n - 5) a'^''fk'^N'^-2n^vua'^ > 0, 

(54) 

so, the larger N, the less anisotropy influences the occur- 
rence of inflation. This is also reflected by the fact that 
f^shear decreases with growing A''. Nevertheless n > 4 is 
still necessary for inflation, no matter how large N is. 

The behaviour of function of n will change 

when N grows. We have seen in the previous section that, 
unlike in the relativistic case, ainf is not a monotonically 
decreasing function of n, and that it has a maximum 
near n = 5. It is not difficult to see that a larger N 
increases the value of Oinf at the maximum and makes 
wider the range of n in which ainf grows with growing 
n. Summarizing, the addition of fields accentuates the 
differences between relativistic and brane cosmologies. 

VI. CONCLUSIONS AND FUTURE 
PROSPECTS 

We have studied Bianchi I type brane cosmologies with 
barotropic perfect fluids, and we have shown that they 
can be interpreted in the language of scalar fields as ex- 
ponential potential models. Prom the expressions here 
obtained, it follows that the models behave, at late times, 
as if the potential had only one exponential term. Consis- 
tently, the differences between the brane and relativistic 
models get washed off as the models evolve in time. The 
effects of the anisotropy disappear in a similar fashion as 
well. 

Our main focus point has been the average evolution 
of the models, in comparison to their relativistic counter- 
parts. In particular, we have studied the features of the 
inflationary behaviour that can be traced back to extra 
dimensional effects but not to the bulk. 



It is a well-known fact that, in the relativistic context, 
the condition 3p < —p (which corresponds to n > 4) is 
necessary and sufficient for inflation to occur. Expressed 
in terms of the potential, it translates into the require- 
ment that the slope of the asymptotically dominating 
potential be smaller than \/2k. In brane cosmology, that 
condition remains necessary for inflation, although we 
show that it is not sufficient. 

Using numeric means wc have found some other new 
interesting results in connection with the kinematics of 
the models. In the relativistic case, the value of the av- 
erage expansion factor at the beginning of inflation de- 
creases monotonically as the ratio p/p becomes more neg- 
ative. In the brane models, the behaviour is basically 
the same, except for the vicinity of 3p < —2/3, where the 
monotonic decrease turns into monotonic increase, until 
a local maximum in the function is reached. Moreover, 
the value of the expansion rate at that extreme point 
does not depend on the amount of anisotropy. 

Taking into account the tight restrictions to the pa- 
rameters of the models set by physical arguments, our 
calculations show that the effect of the anisotropy in the 
evolution will only be noticeable if inflation begins very 
late, which is not a desirable situation. 

We have discussed a possible realization of the assisted 
inflation proposal, which is miderepresented in the liter- 
ature on brane cosmology. In particular, we have consid- 
ered generalizations to multi-field configurations of non- 
interacting scalar fields. The conclusion reached is that 
the aforementioned requirement on the slope of the po- 
tential still applies. On the other hand, an increase in the 
number of fields will make the Hubble factor grow, so 
that inflation will be likelier This means that the brane 
version of the assisted inflation proposal is also feasible. 
Moreover, as A'' grows the models become more isotropic. 

A possible extension to this work would be exploring 
the possibility of generalizing these solutions to include 
bulk effects. It would also be interesting to carry out 
similar studies for anisotropic spacetimes with non-trivial 
space-curvature. Bianchi V type models would be the 
obvious candidate for a first attempt because they are 
the simplest generalization of the open FRW models. 
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